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This  study  presents  the  integration  of  direct  ammonia  solid  oxide  fuel  cell  with  a  gas  turbine  (DA-SOFC / 
GT)  in  a  novel  combined  cooling,  heating  and  power  (CHCP)  cycle.  The  integration  strategy  is  compared 
for  oxygen  ion-conducting  solid  oxide  fuel  cells  (SOFC-O)  and  hydrogen  proton-conducting  solid  oxide 
fuel  cells  (SOFC-H).  Unlike  hybrid  SOFC-absorption  heat  pumps,  the  current  system  is  designed  to  exploit 
the  refrigeration  properties  of  ammonia  to  provide  cooling  with  minimal  bearing  on  complexity  and 
capital  cost.  A  system  analysis  is  developed  to  cover  both  electrochemical  and  thermodynamic  model¬ 
ling.  A  detailed  parametric  study  is  also  conducted  to  investigate  the  effects  of  varying  the  operating 
conditions  and  parameters  on  the  energy  and  exergy  efficiencies  and  the  overall  system  performance. 
The  results  reveal  that  the  SOFC-H  integrated  system  offers  better  performance  than  that  with  the  SOFC- 
O  option.  At  an  operating  temperature  of  1073  K  and  a  pressure  of  500  kPa,  the  respective  energy  and 
exergy  efficiencies  of  the  SOFC-H  integrated  cycle  reach  81.1%  and  74.3%  compared  to  76.7%  and  69.9%  for 
the  SOFC-O.  Under  the  same  operating  conditions,  the  cooling  effectiveness  of  the  system  with  both  SOFC 
types  becomes  14.7%  based  on  the  lower  heating  value  (LHV)  of  ammonia. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Despite  the  economic  recession,  the  global  market  of  SOFC 
systems  is  experiencing  a  robust  growth  with  the  expectation  to 
increase  from  US$360  million  in  2010  to  US$530  million  in  2016  [1  ]. 
This  can  be  partly  attributed  to  the  growing  awareness  of  the 
demonstrated  high  efficiency  and  the  proven  reliability  of  SOFC 
systems  that  position  them  as  ideal  candidates  for  applications  of 
small  to  medium  scale  power  generation  and  public  transportation. 
The  market  opportunity  can  be  further  expanded  by  overcoming 
key  technology  barriers  including  cost  of  materials,  fuel  efficiency 
and  effective  hydrogen  storage. 

In  our  previous  study  [2],  we  examined  the  viability  of  using 
ammonia  as  a  hydrogen  carrier.  The  study  also  included  a  para¬ 
metric  analysis  to  compare  the  performance  of  two  types  of  direct 
ammonia  solid  oxide  fuel  cells  (DA-SOFC)  based  on  the  electrolyte 
used  (i.e.  SOFC-O  or  SOFC-H).  The  findings  showed  that  both 
types  of  the  DA-SOFC  performed  very  well  although  hydrogen 
proton-conducting  DA-SOFC  demonstrated  slightly  better  perfor¬ 
mance  under  all  closed-circuit  conditions.  Accordingly,  it  was 
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concluded  that  ammonia  is  a  highly  efficient  hydrogen  storage 
medium  and  a  suitable  fuel  for  SOFC  applications. 

While  the  scientific  community  and  industry  continue  to 
actively  work  on  improving  the  SOFC  efficiency  and  developing 
new  SOFC  materials,  the  synergy  of  combined  heat  and  power 
(CHP)  systems  plays  a  vital  role  in  the  commercial  viability  and 
economic  sustainability  of  SOFC  technology. 

A  recent  review  conducted  by  Zhang  et  al.  [3]  presented  many 
system  integration  schemes  for  coupling  SOFC  stacks  with  various 
bottoming  cycles  depending  on  the  application  of  use.  One  inter¬ 
esting  application  discussed  in  the  review  is  the  integration  of  SOFC 
with  absorption  refrigeration  system  for  the  combined  delivery  of 
heating,  cooling  and  power  (CHCP).  In  spite  of  the  versatility  of  such 
system,  its  main  disadvantages  include  large  volume,  low  cooling 
efficiency  and  high  capital  cost  due  to  large  number  of  components. 

Zamfirescu  and  Dincer  [4]  proposed  a  hydrogen-fed  integrated 
SOFC/GT  system  for  vehicular  applications.  The  performance  of  the 
system  was  assessed  based  on  the  optimization  objective  of 
maximum  efficiency  or  maximum  power  density.  Additionally,  an 
exergy  analysis  was  carried  out  to  outline  the  major  sources  of 
exergy  destruction  in  the  system.  It  was  found  that  the  fuel  cell 
stack  can  occupy  about  75%  of  the  total  system  volume  and  reduce 
the  compactness  of  the  system  by  40%  when  the  maximum  system 
efficiency  is  achieved.  They  also  concluded  that  the  largest  exergy 
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destruction  occurred  at  the  fuel  cell  stack  and  accounted  for  60%  of 
the  total  exergy  destruction  in  the  system. 

Finally,  another  study  by  Dincer  et  al.  [5]  focused  on  the  exergy 
analysis  of  an  elaborate  SOFC/GT  system  fuelled  with  methane.  The 
study  highlighted  the  advantage  of  higher  operating  voltage  on  the 
thermal  and  exergy  efficiency  of  the  system.  It  was  also  concluded 
that  the  highest  thermodynamic  irreversibilities  emerged  from  the 
electrochemical  and  combustion  processes.  In  particular,  37.1%  of 
the  total  exergy  destruction  occurred  in  the  fuel  cell  stack  while 
34.8%  of  the  total  exergy  was  destroyed  in  the  combustion  chamber. 

The  objective  of  this  work  is  to  propose  integrated  SOFC  systems 
using  ammonia  as  a  direct  fuel.  A  comparative  assessment  of 
oxygen  ion-conducting  solid  oxide  fuel  cells  (SOFC-O)  and  hydrogen 
proton-conducting  solid  oxide  fuel  cells  (SOFC-FI)  is  conducted.  In 
order  to  demonstrate  and  evaluate  these  systems  comprehensive 
energy  and  exergy  analyses  are  conducted,  and  their  energy  and 
exergy  efficiencies  are  assessed  for  practical  applications.  Also, 
energy  and  exergy  efficiencies  of  a  novel  DA-SOFC/GT  system  are 
examined.  Furthermore,  a  detailed  parametric  study  is  carried  out 
to  investigate  the  effects  of  design  and  operating  conditions  on  the 
performance  of  the  system. 

2.  Description  and  operational  principle  of  the  DA-SOFC 

Regardless  of  the  type  of  electrolyte  material,  fuel  in  the  form  of 
gaseous  ammonia  is  subjected  to  catalytic  dehydrogenation  as  it  is 
fed  through  the  anode  channel  of  the  DA-SOFC.  The  extracted 
hydrogen  is  then  electrochemically  oxidized  in  the  fuel  cell  to 
produce  electric  power. 

Depending  on  the  thermodynamic  equilibrium  of  the  ammonia 
decomposition  reaction  at  a  given  temperature  and  pressure, 
nitrogen  gas  may  be  produced  in  high  concentration.  The  inert 
nitrogen  can  dilute  the  hydrogen  concentration  at  the  anode 
resulting  in  reduction  of  the  reversible  cell  potential.  The  nitrogen 
produced  is  discharged  from  the  anode  compartment  along  with 
other  gaseous  effluents. 

Fig.  1  depicts  the  configuration  of  the  DA-SOFC  with  oxygen 
ion-conducting  electrolyte.  Oxygen  from  air  is  transported  across 
the  cathode  layer  to  the  cathode-electrolyte  interface  where  it  is 
reduced  to  oxygen  ion.  The  oxygen  ions  are  then  transported 
through  the  electrolyte  layer  to  electrochemically  react  with 
hydrogen  at  the  anode-electrolyte  interface.  The  water  vapour 
produced  by  the  reaction  migrates  to  the  anode  channel  thereby 
further  reducing  the  hydrogen  concentration.  The  water  vapour 
exits  the  fuel  cell  along  with  unreacted  fuel  and  other  effluents. 

The  arrangement  of  the  DA-SOFC  with  hydrogen 
proton-conducting  electrolyte  is  shown  in  Fig.  2.  Hydrogen  is 
oxidized  at  the  anode-electrolyte  interface  into  protons  and  elec¬ 
trons.  The  protons  are  transported  through  the  electrolyte  layer  to 
the  cathode-electrolyte  interface  to  undergo  electrochemical 
reaction  with  oxygen  from  the  cathode  side.  The  water  vapour 
produced  by  the  fuel  cell  reaction  is  then  transferred  to  the  bulk  gas 
in  the  cathode  channel  where  it  exits  the  cell  with  other  effluents. 


NH3 


NH3  3/2  H2  ♦  1/2  N2 


=>  H2.NiNH3.  H20 


Anode 

Electrolyte 

Cathode 


h2o  \\2o 

♦h2 

To*  |o2  o*T 

ts, _ Oil 

00 


02,  n2 


=>  02.N2 


NHj  t=>  NHj  ->  3/2  H 2  +  1/2  Nj 


=>  H2.N2.NII, 


Anode 

Electrolyte 

Cathode* 


h2*  »h2 

1 

h4 

£ 

°2  1 

\  |h2o  1 

k°2 

■""H 

1 - - - 1 

V 


02,N2  «=> 


o2.n2,h2o 


Fig.  2.  Schematic  representation  of  hydrogen  proton-conducting  DA-SOFC. 


This  is  a  major  advantage  since  hydrogen  in  the  anode  compart¬ 
ment  is  not  diluted  by  the  water  vapour  generated  from  the  fuel  cell 
reaction. 

3.  Description  of  integrated  DA-SOFC/GT  system 

The  DA-SOFC  can  be  integrated  with  a  bottoming  thermody¬ 
namic  cycle  where  the  fuel  cell  is  incorporated  into  a  gas  power 
system  in  a  CHCP  system.  In  the  particular  arrangement  depicted  in 
Fig.  3,  the  fuel  tank  contains  pressurized  ammonia  for  on-demand 
delivery  to  the  system  as  needed.  As  the  ammonia  is  extracted 
during  the  operation  of  the  system,  its  temperature  and  pressure 
decrease  by  an  amount  proportional  to  its  enthalpy  of  vaporization. 
In  order  to  maintain  the  tank  pressure  at  the  desired  level, 
a  closed-loop  heat  exchanger  (coil)  is  utilized  to  transfer  rejected 
heat  from  a  suitable  fluid  to  the  ammonia  in  the  tank.  As  an  added 
benefit,  the  heat  transfer  process  and  the  refrigeration  properties 
of  ammonia  are  exploited  to  create  a  cooling  effect  in  the  medium 
flowing  inside  the  closed-loop  heat  exchanger.  This  configuration 
provides  a  simple  and  cost-effective  alternative  to  combined 
SOFC-Absorption  refrigeration  cycles. 

The  extracted  ammonia  is  further  heated  in  a  two-stage  heat 
exchanger  arrangement  (HX1  and  HX2)  to  achieve  a  reasonable 
temperature  before  entering  the  fuel  cell  stacks.  At  this  point,  air  is 
drawn  into  the  compressor  unit  to  be  compressed  to  the  desired 
pressure.  Although  the  temperature  of  air  rises  during  the 
compression  process,  further  heating  is  achieved  by  passing  the 
compressed  air  through  a  heat  exchanger  (HX3)  upstream  of  the 
fuel  cell  stacks.  The  fuel  and  oxygen  from  air  (oxidant)  are  elec¬ 
trochemically  reacted  inside  the  fuel  cells  to  produce  electric  power 
while  the  heat  generated  due  to  thermodynamic  irreversibilities  is 
used  for  the  thermal  decomposition  of  ammonia  and  additional 
heating  of  the  reactants  and  products. 

As  governed  by  the  fuel  utilization  factor,  a  fraction  of  the 
hydrogen  dissociated  from  the  ammonia  exits  the  fuel  cell  stacks 
along  with  other  gaseous  effluents  to  be  burned  in  the  adiabatic 
combustion  chamber.  The  temperature  of  the  gaseous  mixture 
leaving  the  combustion  chamber  is  sufficiently  high  for  expansion 
and  production  of  useful  mechanical  power  in  a  gas  turbine.  Some 
of  the  power  produced  is  used  to  drive  the  compressor  while  the 
remainder  is  used  to  drive  an  electric  generator.  The  gas  mixture 
leaving  the  turbine  is  split  into  two  separate  branches  that  feed  the 
heat  exchangers  discussed  earlier.  Finally,  heat  is  recovered  from 
the  exhaust  gas  which  leaves  the  system  through  the  recuperator  at 
a  moderate  temperature. 


4.  Thermodynamic  analysis 

4.1.  Thermodynamic  equilibrium  of  ammonia  decomposition 


When  heated,  ammonia  decomposes  to  hydrogen  and  nitrogen 
Fig.  1.  Schematic  representation  of  the  oxygen  ion-conducting  DA-SOFC.  according  to  the  following  reaction: 
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Fig.  3.  Schematic  diagram  of  integrated  DA-SOFC/GT  system. 


nh3-|h2  +  1n2  (i) 

The  above  reaction  is  endothermic  and  depends  on  the 
temperature  and  pressure  of  the  system  among  other  factors.  The 
reaction  proceeds  to  the  right  as  the  temperature  increases  or  to 
the  left  as  the  pressure  increases  thus  producing  less  hydrogen.  The 
thermodynamic  equilibrium  of  ammonia  can  be  obtained  by 
minimizing  the  total  Gibbs  energy  of  the  reaction: 

(AGsystem)^  p  =  0  (2) 


Csystem  -  (^>  [gg  +  Rrin(y,P)]  ^  +  (£ n4)condensed  (5) 

As  an  advantage  of  this  method,  chemical  reactions  need  not  to 
be  specified.  Instead,  Gibbs  energy  minimization  (GEM)  requires 
data  of  the  standard  Gibbs  function  of  formation  of  all  chemical 
species  which  are  expected  to  be  involved  or  contribute  to  the 
thermodynamic  equilibrium  under  consideration. 

4.2.  Cell-level  thermodynamic  and  electrochemical  modelling  of 
SOFC 


where  the  Gibbs  energy  of  the  system  can  be  defined  as  the 
summation  of  the  product  of  moles  of  chemical  species  i  multiplied 
by  the  respective  specific  Gibbs  energy  at  constant  temperature  and 
pressure 

^system  =  njgj  (3) 

The  specific  Gibbs  energy  can  be  described  as  the  sum  of  the 
standard  Gibbs  function  of  formation  and  the  chemical  activity  of 
a  given  species  i.  For  real  gases,  the  activity  of  species  is  defined  as 
the  ratio  of  its  fugacity  in  the  system  to  that  at  standard  tempera¬ 
ture  and  pressure  (STP)  conditions: 

I,  =  gfi  +  RT  In  j‘Q  (4) 

h 

If  ideal  or  perfect  gas  conditions  are  assumed,  the  activity  of  the 
gaseous  species  is  equal  to  its  partial  pressure  in  the  system.  For 
condensed  matter  (solid  or  liquid),  the  activity  is  equal  to  unity. 
Therefore,  the  total  Gibbs  energy  of  the  system  becomes 


The  equilibrium  potential  of  any  fuel  cell  is  affected  by  the 
partial  pressures  of  the  chemical  species  involved  in  the  electro¬ 
chemical  reaction.  Under  open-circuit  conditions,  the  potential 
(voltage)  can  be  written  as 


E  = 


fo  +  gm 


Ph2(Po2)1/2 

Ph2o 


(6) 


where 


E° 


-AG 


(7) 


Furthermore,  the  examination  of  the  potential  of  solid  oxide  fuel 
cells  at  closed-circuit  conditions  requires  the  characterization  of 
the  over-potentials  affecting  the  cell  operation.  The  working 
potential  or  voltage  of  SOFC  can  be  characterized  as 


V  —  E  <Pact,an  ^act,ca  VQ.  ^conc,an  ^conc,ca 


(8) 
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The  activation  over-potential  is  an  irreversible  loss  related  to  the 
electrode  kinetics  and  represents  the  difference  between  the  actual 
potential  (closed-circuit  conditions)  and  the  reversible  potential 
(open-circuit  conditions).  The  implicit  relationship  between  the 
current  density  and  the  activation  over-potential  can  be  described 
using  the  Butler— Volmer  expression  [6] 


Df  =  0.0018583- 


r3|E  +  _l 


'\  MjJ  Pa}flDM 


(13) 


and  the  mean  characteristic  length  (cry)  of  the  molecular  collision 
diameter  of  species  i  and  j  is  given  by 
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The  exchange  current  density  (Jo)  is  a  measure  of  the  magnitude 
of  electron  activity  at  the  equilibrium  potential  of  the  electrode.  Its 
value  is  highly  influenced  by  the  electrode  material,  structure  and 
other  factors  such  as  the  temperature  of  the  reaction  and  the  length 
of  the  triple-phase  boundary  (TPB)  [7].  The  charge  transfer  coeffi¬ 
cient  (a)  describes  the  effect  of  the  electric  potential  on  the  ratio  of 
the  forward  to  the  reverse  activation  barrier.  For  most  electro¬ 
chemical  reactions  in  solid  oxide  fuel  cells,  this  value  is  assumed  to 
be  0.5  [8].  As  such,  the  explicit  relationship  of  the  activation  over¬ 
potential  is 


RT  (  1  \ 

<Pact,y  =  ~Ys^n  ft-1  izj—)  ’  7  =  anode,  cathode  (9b) 

The  electronic  resistance  of  the  fuel  cell  electrode  and  inter¬ 
connect  is  usually  very  small  and  is  considered  negligible  when 
compared  to  the  ionic  resistance  of  the  electrolyte.  The  ionic 
resistance  can  be  described  by  Ohm’s  law  [9]  as 

n=mQ  (io) 

Here,  the  electrolyte  resistance  (Pq)  is  a  function  of  the  electrolyte 
properties  and  is  highly  dependent  on  the  operating  temperature. 

The  concentration  over-potential  accounts  for  losses  incurred  by 
the  resistance  of  the  porous  electrode  to  the  transport  of  gaseous 
species  between  the  gas  channel  and  the  reaction  sites  at  TPB.  The 
mass  transport  in  the  electrodes  is  driven  by  the  diffusion  of 
reacting  species  due  to  concentration  gradient  as  well  as  the 
permeation  caused  by  the  pressure  gradient.  Several  mass  trans¬ 
port  models  with  varied  accuracies  have  been  used  to  capture  the 
effects  of  concentration  over-potential  on  the  performance  of 
fuel  cells. 

It  has  been  previously  demonstrated  that  the  highest  accuracy 
can  be  achieved  using  the  dusty  gas  model  (DGM)  [10]  and  the 
mean  transport  pore  model  (MTPM)  [11].  In  particular,  DGM  visu¬ 
alizes  the  porous  medium  as  a  collection  of  large  spherical  particles 
suspended  in  space  by  external  forces.  The  model  takes  into 
account  the  effects  of  diffusion  and  permeation  mechanisms  on  the 
concentration  of  reacting  species  which  is  essential  when  model¬ 
ling  pressurized  fuel  cell  systems.  The  one-dimensional  multi- 
component  mass  transport  DGM  can  be  written  as 
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where  the  effective  Knudsen  diffusion  accounts  for  the  porosity  and 
tortuosity  of  the  electrode  such  that  [12]: 


(14) 


The  collision  integral  (Qd,  y)  is  a  function  of  temperature  and  the 
Lennard -Jones  parameter  (k/ey)  such  that 


n  _  1.06036  0.19300  1.03587 

DM  -  p0. 15610  +  exp(0.47635r)  +  exp(1.52996r) 
1.76474 

+  exp(3.8941ir) 

where 


r  =  —  (16) 

£ij 

and 

£ij  =  JziZj  (17) 

The  values  of  the  parameter  in  Eqs.  (12)— (17)  have  been 
obtained  from  references  [13]  and  [14].  The  permeability  of  the 
porous  electrode  is  estimated  using  the  Kozeny-Carman  relation¬ 
ship  [15] 


Bo 
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Furthermore,  the  viscosity  of  gas  mixture  can  be  calculated 
using  the  following  semi-empirical  equation  [16]: 


„  _  TJiPiVMi 

E  yiy/Mi 


(19) 


The  gas  viscosities  are  obtained  from  reference  [13].  As  pre¬ 
sented  by  Zhu  et  al.  [17],  the  summation  of  Eq.  (11)  enables  the 
evaluation  of  the  pressure  gradient  (dP/dx)  across  the  electrode 
layer  as 


At  TPB,  the  molar  flux  of  the  gaseous  reactants  involved  in  the 
electrochemical  reaction,  namely  hydrogen  and  oxygen,  can  be 
related  to  the  current  density  as 


Ni  =  ^  (2D 

Depending  on  the  type  of  electrolyte  (SOFC-O  or  SOFC-H),  the 
molar  flux  of  water  vapour  can  be  calculated  using  Graham’s  law  of 
diffusion  which  governs  the  diffusion  process  of  gas  mixtures  [10]: 
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The  effective  binary  diffusion  of  chemical  species  can  be 
obtained  using  the  Chapman-Enskog  equation  [13]: 


EN<VMi  =  0  (22) 

i 

The  molar  fluxes  of  all  other  non-reacting  species  are  equal  to 
zero  [18].  Finally,  the  relationship  between  the  concentration 
over-potential  and  the  partial  pressures  can  be  written  as  [19] 
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(23c) 

(23d) 


\PH2ofpToP2BJ 

4.3.  Fuel  and  oxidant  utilization 

The  utilization  of  ammonia  can  be  defined  in  terms  of  the  actual 
supply  and  consumption  of  the  fuel  or  its  hydrogen  equivalent  such  that 


Uf  = 


(Fuel)consumed  _  (H2)consumed 
(Fuel)supplied  (^2)  supplied 


(24) 


Similarly,  the  utilization  of  air  can  be  expressed  as  a  function  of 
the  air  supply  and  usage  or  its  oxygen  (oxidant)  equivalent  since  air 
can  be  assumed  to  be  composed  of  79%  nitrogen  and  21%  oxygen  by 
volume. 


Uo  = 


(Air) consUmed  _  (^2) consumed 


(Air)  supplied  (O2  Supplied 
4.4.  Energy  and  exergy  analyses  of  integrated  system 


(25) 


The  exergy  destruction  due  to  thermodynamic  irreversibilities  in 
any  component  can  be  written  as 


F*D,cmP  —  Qj 


l-4£ 


-W  +  ^(EXi)  in -£(**, 


i/OUt 


1  1 
and  the  total  exergy  destruction  in  the  system  is 
F^d, total  —  I>D,anp 


(31) 


(32) 


The  electric  energy  efficiency  is  dependent  on  the  lower  heating 
value  of  the  fuel 


Ve  = 


W t 


AC  ' 


Wr 


h  LHV 


and  the  electric  exergy  efficiency  can  be  written  as 


ie  = 


hex°ch 


The  total  energy  efficiency  of  the  system  is 

=  Wac  +  WG  +  Qrec 
Vt  ~  hLHV  +  Qcl 

and  the  total  exergy  efficiency  of  the  system  is 
WAc  +  WG  + 


it  = 


nex°ch  +  EXq£L 


(33) 


(34) 


(35) 


(36) 


4.5.  Energy  and  exergy  analyses  of  system  components 

In  this  section,  ideal  gas  conditions  are  assumed.  Unless  other¬ 
wise  stated,  all  components  and  devices  are  considered  to  be  well 
insulated  to  minimize  heat  losses  to  negligible  levels.  Furthermore, 
pressure  losses  have  been  ignored  by  virtue  of  component  design. 

•  Fuel  cell  stack 


At  steady  state  operation,  the  energy  balance  of  the  system  or 
any  of  its  components  can  be  written  as 

Q-  -  W  +  £  (fiihi) in  -  £  (nfii) out  =  0  (26) 

i  i 

The  exergy  associated  with  material  flow  is  the  sum  of  the 
physical  and  the  chemical  exergy 


The  thermodynamic  processes  taking  place  inside  the  fuel  cell 
module  are  depicted  in  Fig.  4.  The  ratio  of  hydrogen  to  ammonia 
entering  the  fuel  cell  is  determined  by  the  thermodynamic  equi¬ 
librium  and  depends  on  the  temperature  and  pressure  of  the  cor¬ 
responding  reaction.  Therefore,  the  heat  generated  by  the 
isothermal  electrochemical  reaction  of  hydrogen  and  oxygen  inside 
the  fuel  cell  based  on  the  molar  flow  of  ammonia  is 


Ex  =  Ex  ph  +  Exch 

(27) 

where 

EXph  =  TLhi  [fa  -  So)  -  To  (Si  -  So)] 
i 

(28) 

Exch  =  n  +  lny,- 

(29) 

i  i 


The  previous  expressions  assume  negligible  change  in  kinetic 
and  potential  energies.  In  addition,  the  exergy  of  work  is  work 
itself;  however,  the  exergy  delivered  by  heat  transfer  is  governed  by 
the  Carnot  efficiency  such  that 


Fxq  =  Qj 


Jo 

T 


(30) 


In  this  study,  the  exergy  is  evaluated  against  reference 
temperature  and  pressure  of  298  K  and  101.325  kPa,  respectively. 


Q  —  xflUfTp  gAs  +  I  (<Pact,an  +  ^act,ca  +  ^ohm  +  cone, an  +  ^conc,an) 

(37) 

The  generated  heat  is  not  only  used  for  the  internal  thermolysis 
of  fuel  but  also  to  supply  sensible  heat  to  the  reactants  and  prod¬ 
ucts.  Furthermore,  the  temperature  of  all  effluents  is  assumed 
uniform  and  equal  to  the  fuel  cell  temperature.  All  surplus  heat  is 
rejected  outside  the  fuel  cell  stacks.  The  thermal  efficiency  of  the 
fuel  cell  stack  can  be  written  as 


^FC  = 


We 

A H 


and  the  exergy  efficiency  is 
F*d,fc 


foe  1- 


(38) 


(39) 


Zifaexi) in  -  E ifaexOout 
Finally,  the  actual  AC  power  produced  by  the  fuel  cell  stack  is 
^ac  =  ^invWfc  (40) 
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V'T 


WT 

£;(' i;exi)in  -  £;(' iiexOout 


(48) 


•  Air  compressor 

This  component  receives  air  at  ambient  temperature  and  pres¬ 
sure  and  compresses  it  to  the  desired  pressure.  This  also  results  in 
an  increase  in  the  air  temperature  at  the  outlet  of  the  compressor: 

Tout  =  Tin(PR)(,<-1)/k  (49) 


Fig.  4.  Thermodynamic  processes  of  the  fuel  cell  stack. 


where 


•  Combustion  chamber 

Unreacted  hydrogen  and  other  effluents  including  excess  air 
leave  the  fuel  cell  stacks  and  enter  the  combustion  chamber  where 
no  work  is  done.  It  is  also  assumed  that  excess  air  ensures  complete 
combustion.  Applying  the  principles  of  energy  and  mass  conser¬ 
vation  on  the  component,  yields 

E  ("lSl)in  =  E  ("AOout  (41) 

i  i 

For  the  current  system,  this  translates  to 

aH2  +  b02  +  dH20  +  kN2  -►  (a  +  d)H20  +  (b  -  02  +  kN2  (42) 

Knowing  the  molar  flow  rate  and  the  chemical  composition  of 
the  combustion  reactants  and  products,  the  adiabatic  combustion 
temperature  can  be  iteratively  determined.  Finally,  the  exergy 
efficiency  of  the  combustion  chamber  is 


(50) 


ZiyiCvj 


(51) 


Similar  to  the  gas  turbine,  the  enthalpy  of  air  at  isentropic 
conditions  can  be  related  to  its  actual  enthalpy  through  the 
isentropic  efficiency  of  the  compressor  such  that 


Ei(rciftjs,out  -  £ ffl)jn 

J2i(^Mout  —  X)i(d(dj)in 


from  which  the  actual  outlet  temperature  of  air  is  iteratively 
obtained.  The  exergy  efficiency  of  the  compressor  is 

=  E,(niex,)out  -  £,(n,ex,)in  (53) 

wc 


fc:  1 


E*d,cc 

£<(£*)  in 


•  Gas  turbine 


•  Electric  generator 

The  surplus  power  representing  the  difference  between  the 
mechanical  power  generated  by  the  turbine  and  that  consumed  by 
compressor  can  be  captured  using  the  back  work  ratio  (BWR) 


When  the  hot  gaseous  mixture  leaves  the  combustion  chamber 
and  enters  the  gas  turbine,  it  expands  to  deliver  useful  mechanical 
power.  With  the  knowledge  of  the  inlet  temperature,  the  outlet 


temperature  can  be  determined  by 

Tout  =  Tin(  pjj) 

(44) 

where 

P 

PR  = 

Pout 

(45) 

k  _  TsiyiCpj 

£iy.-4,i 

(46) 

The  actual  outlet  temperature  of  the  turbine  is  higher  due  to  the 
effects  of  thermodynamic  irreversibilities.  The  isentropic  efficiency 
of  the  turbine  can  be  applied  to  account  for  the  difference  such  that 


Wr 

BWR  =  (54) 

Wj 

It  is  definitely  desired  to  maintain  this  ratio  as  low  as  possible. 
The  surplus  power  is  used  to  drive  an  electric  generator: 

=  Vg(Wt-Wc)  (55) 

•  Fleat  exchangers 

Multiple  gas-to-gas  heat  exchangers  are  utilized  in  this  system. 
The  effectiveness  method  [20]  is  used  to  determine  the  outlet 
temperature  of  the  respective  gaseous  stream  such  that 

Q.  =  ?^C,min  (^h,in  —  ^c,in)  (56) 

where  Hc<  min  is  the  minimum  heat  capacity  rate  of  the  hot  or  the 
cold  stream.  As  such,  the  outlet  temperature  of  the  hot  stream  is 
given  as 


v  =  SiQVhQ  in  SiQVhQout 

S,T  £K"Mn-£i(^i)s,out 

The  actual  temperature  of  gas  mixture  leaving  the  turbine  can 
be  iteratively  determined  from  the  actual  enthalpy  obtained.  The 
exergy  efficiency  of  the  turbine  can  be  written  as 


Q  —  (^icP,i)h  (Thfn  ^h, out) 

i 

and  the  cold  stream 

Q.  =  (^P,i)c(Tc,°ut  —  ^c,in) 

i 


(57) 


(58) 
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•  Mixing  chamber 


Gas  streams  of  different  molar  flow  rates  and  temperatures 
flowing  through  different  branches  are  mixed  in  a  mixing  chamber 
before  entering  heat  exchanger  (HX1).  Under  the  assumption  of 
adiabatic  mixing  of  gases,  the  temperature  of  the  gaseous  mixture 
exiting  the  mixing  chamber  can  be  obtained  from 

=  EW)out  (59) 

i  i 

The  exergy  destruction  can  be  determined  from  Eq.  (31)  while 
the  exergy  efficiency  can  be  written  as 


^MC  =  1 


£*d,mc 

£i(ni®*i)in 


(60) 


5.  Results  and  discussion 


Fig.  5.  Thermodynamic  equilibrium  of  ammonia. 


The  minimization  of  Eq.  (5)  was  performed  using  the  optimi¬ 
zation  function  in  MATLAB  (Math Works,  Natick,  MA,  USA)  [21].  In 
order  to  obtain  the  molar  fractions  and  partial  pressures  of  the 
reacting  species  at  TPB,  a  computer  code  was  developed  using 
function  (ode45)  in  MATLAB  to  numerically  solve  the  set  of 
simultaneous  differential  equations  represented  by  Eq.  (11). 
Furthermore,  the  required  thermodynamic  and  chemical  exergy 
data  were  obtained  from  references  [22-24].  All  thermodynamic 
calculations  have  been  performed  using  EES  (F-Chart  Software, 
Madison,  WI,  USA)  [25]. 

5  A.  Thermolysis  of  ammonia  in  the  DA-SOFC 


electrodes  and  the  electrolyte  as  well  as  the  microstructural 
properties  of  the  cell  components  which  are  highly  dependent  on 
the  processing  method  as  discussed  by  Ringuede  et  al.  [30].  The 
results  in  Fig.  6  show  excellent  agreement  (within  2%)  between  the 
numerical  simulation  and  the  experimental  data. 

It  is  also  worth  noting  that  experimental  tests  are  typically 
performed  with  low  utilization  of  reactants  in  order  to  avoid  the 
risk  of  local  overheating  and  material  degradation  which  can  occur 
at  high  current  densities.  In  practice,  higher  utilization  factors  are 
always  desirable  despite  the  lower  efficiency  and  specific  power 
output  of  the  fuel  cell  stack. 


The  use  of  the  GEM  method  revealed  that  the  decomposition  of 
ammonia  is  favoured  at  high  temperature  and  low  pressure  as 
shown  in  Fig.  5.  In  particular,  full  conversion  of  ammonia  can  be 
achieved  at  a  temperature  of  700  I<  and  atmospheric  pressure. 
Nonetheless,  it  is  well  known  that  the  kinetic  rate  of  reaction  is 
rather  slow  and  often  requires  the  use  of  a  suitable  catalyst  to 
promote  efficient  conversion  of  ammonia.  Complete  conversion 
may  be  achieved  over  expensive  metal  catalysts  such  as  ruthenium 
and  platinum  at  temperatures  around  600-650  K  [26];  however, 
more  economic  alternatives  include  iron  and  nickel-based  cata¬ 
lysts.  Since  the  operating  temperature  of  solid  oxide  fuel  cells 
commonly  ranges  between  900  and  1400  K  [27],  it  is  safe  to  neglect 
the  kinetic  rate  of  reaction  and  assume  that  full  decomposition  of 
ammonia  has  been  attained  within  the  porous  anode  layer  of  the 
fuel  cell. 

5.2.  Model  validation 

In  the  open  literature,  experimental  investigations  related  to  the 
DA-SOFC  are  scarce  and  rarely  contain  all  the  necessary  information 
and  parameters  needed  for  model  evaluation.  Therefore,  certain 
assumptions  were  made  and  some  model  parameters  were 
adjusted  to  agree  with  the  experimental  data  as  reflected  in  Table  1. 
Ma  et  al.  [28]  conducted  an  experiment  in  which  ammonia  was 
used  as  a  fuel  in  anode-supported  SOFC-O  with  the  arrangement  of 
Ni-YSZ|YSZ|LSM-YSZ.  The  cell  was  tested  at  atmospheric  pressure 
and  different  temperatures  between  923  K  and  1123  K.  The 
performance  of  the  cell  was  characterized  by  the  high  interfacial 
resistance  which  accounted  for  about  50%  of  the  total  resistance  of 
the  cell.  The  resistance  of  the  electrolyte  was  also  higher  than 
expected  at  the  same  temperature  as  reported  elsewhere  [19,29].  In 
part,  this  can  be  attributed  to  the  poor  contact  between  the 


5.3.  Energy  and  exergy  performance  of  integrated 
DA-SOFC/GT  system 

As  noted  earlier,  it  is  desired  to  operate  the  fuel  cell  at  the 
highest  possible  utilization  factor  while  maintaining  reasonable 
concentration  of  reactants  and  performance  of  the  fuel  cell.  This 
value  may  vary  between  80  and  95%  for  fuels  and  20-35%  for 
oxidants.  Table  2  reflects  the  parameters  used  for  the  numerical 
simulation  of  the  base  case  study.  The  polarization  curve  of 
the  DA-SOFC  is  shown  in  Fig.  7.  The  peak  power  densities  are 
1875  W  m-2  at  0.67  V  for  SOFC-O  and  4507  W  m“2  at  0.58  V  for 


Table  1 

Parameters  used  in  the  validation  of  oxygen  ion-conducting  DA-SOFC. 


Parameter 

Value 

Operating  temperature,  T  (K) 

1123 

Operating  pressure,  P  (kPa) 

101 

Fuel  utilization,  Uf{%) 

5 

Oxidant  (air)  utilization,  U0  (%) 

5 

:wcsf)«p(-ioD 

Exchange  current  density  of 
anode  Jan  (A  m-2) 

7  x  109( 

Exchange  current  density  of  cathode, 

7  x  109( 

?r-(J3r) 

Jca  (A  ITT2) 

Electrode  porosity,  e 

0.4 

Electrode  tortuosity,  t 

4.25 

Average  pore  radius  of  electrode, 

0.5 

rp  (pm) 

Resistivity  of  electrolyte,  Rq  (O  m) 

0.735 

Anode  thickness,  (pm) 

500 

Electrolyte  thickness,  (pm) 

30 

Cathode  thickness,  (pm) 

10 

Sources  [28,31-33]. 
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SOFC-H.  The  significant  difference  of  144%,  based  on  peak  power 
density,  was  mainly  due  to  the  higher  concentration  of  hydrogen  at 
the  anode  of  SOFC-H  compared  to  that  at  the  anode  SOFC-O.  This 
was  also  reflected  by  the  difference  in  the  open  cell  voltage  as  well 
as  the  detrimental  impact  of  the  concentration  polarization  which 
is  most  pronounced  in  the  case  of  SOFC-O. 

Furthermore,  the  respective  fuel  cell  could  not  be  operated  at 
the  corresponding  peak  power  density  due  to  the  overall  arrange¬ 
ment  of  integrated  system  and  the  heating  requirements  of  the 
thermolysis  of  ammonia.  The  reduced  operational  voltages  and 
current  densities  are  listed  in  Table  3. 

The  temperature  distribution  of  the  system  is  depicted  in  Fig.  8 
and  the  state  points  correspond  to  those  in  Fig.  3.  For  the  base  case 
study,  the  temperature  of  the  fuel  cell  effluents  was  1073  I<  (T5  and 
T6)  and  increased  to  1193  I<  (T7)  in  the  combustion  chamber  due  to 
the  adiabatic  combustion  of  unused  fuel.  The  hot  exhaust  gases 
entered  the  turbine  and  expanded  to  produce  mechanical  work  and 
exited  the  turbine  at  853  K  (T8).  At  this  point,  the  exhaust  splits  into 
two  branches  to  deliver  heating  gas  streams  to  the  fuel  and  air  heat 
exchangers. 

The  difference  of  130  K  between  the  temperature  of  the  anode 
(T3)  and  the  cathode  (T4)  inlets  may  potentially  cause  mechanical 
stresses  which  can  lead  to  problems  with  sealing,  leakage  and 
cracking  of  ceramic  components.  This  temperature  mismatch  can 
be  explained  by  realizing  that  95%  of  the  turbine  exhaust  volume 


Table  2 

Parameters  used  for  numerical  simulation  of  the  base  case  study. 


Parameter  Value 

Fuel  utilization  factor,  Uf  [%]  85 

Oxidant  utilization  factor,  U0  [%]  25 

Stack  temperature,  TFC  [I<]  1073 

Stack  Pressure,  PFc  [  kPa ]  505 

Pressure  Ratio,  PR  5 

Cumulative  cell  area,  AFC  [m2]  250 

Compressor  efficiency,  %  c  [%]  85 

Turbine  efficiency,  T  [%]  90 

Heat  exchanger  effectiveness,  ?  [%]  90 

AC  generator  efficiency,  r]G  [%]  95 

DC/A C  inverter,  T]my  [%]  98 

Ambient  temperature,  Tamb  [K]  298 

Ambient  pressure,  Pamb  [kPa]  101.3 

Composition  of  oxygen  in  air,  [%]  21 

Composition  of  nitrogen  in  air,  [%]  79 


Fig.  7.  Polarization  curve  of  oxygen  ion  and  hydrogen  proton-conducting  DA-SOFC 
(base  case). 

was  diverted  as  shown  in  Fig.  3  to  heat  exchanger  (HX3)  through 
branch  (8a)  in  order  to  satisfy  the  high  sensible  heat  demands  of 
the  incoming  air.  It  also  resulted  in  a  low  heat  transfer  rate  in  (HX2) 
that  could  not  raise  the  temperature  of  the  fuel  to  match  that  of  the 
air  entering  the  fuel  cell. 

The  total  specific  electric  power  generated  by  the  system  with 
SOFC-H  option  was  about  8%  higher  than  that  of  the  system  with 
SOFC-O  as  given  in  Table  4.  The  difference  corresponds  to  the  lower 
exergy  destruction  in  the  system  with  SOFC-H  option  due  to 
reduced  thermodynamic  irreversibilities.  The  tabulated  data  also 
revealed  that  more  than  62%  of  the  total  power  produced  is  deliv¬ 
ered  by  the  fuel  cell  stacks.  This  reasonable  performance  can 
be  further  improved  by  reducing  the  exergy  destruction  in  the 
fuel  cell. 

In  order  to  maintain  the  pressure  of  the  ammonia  exiting  the 
tank  at  505  kPa  (with  a  saturation  temperature  of  277.3  K),  heat  was 
transferred  from  the  closed-loop  coil  to  the  ammonia  in  the  tank  at 
a  rate  of  46.7  kj  mol-1  of  ammonia.  The  heat  rejection  produced 
a  cooling  effect  in  the  fluid  passing  through  the  coil  with  an 
effectiveness  of  14.7%  based  on  the  LHV  of  ammonia. 

Fig.  9  shows  the  exergy  destruction  within  a  system  component 
as  a  percentage  of  the  total  exergy  destruction  of  the  power  plant.  It 
can  be  seen  that  the  distribution  of  exergy  destruction  is  very 
similar  in  both  systems.  The  highest  exergy  destruction  took  place 
in  the  fuel  cell  stacks  due  to  thermodynamic  irreversibilities  related 
to  internal  electrochemical  processes  as  well  as  external  losses  due 
to  the  large  thermal  gradient  between  the  fuel  cell  stack  and  its 
surroundings.  Thus,  the  figure  reiterates  that  efforts  should  be 
focused  on  high  exergy  destruction  sites  like  the  fuel  cell  module  to 
minimize  its  thermodynamic  losses  and  achieve  maximum 
performance  and  efficiency. 


Table  3 

Operational  voltages  and  current  densities  of  the  integrated  system  at  different 
conditions. 


Fuel  cell 

SOFC-O 

SOFC-H 

Parameters 

V[V] 

7  [A  m  2] 

V[V] 

J  [Am  2] 

Base  case 

0.52 

3185 

0.59 

7640 

TFC  =  1173  I< 

0.44 

3545 

0.52 

7935 

TVc  =  1273  I< 

0.36 

3440 

0.45 

8200 

Uf  =  0.90 

0.48 

2190 

0.55 

6325 

Uf=  0.95 

0.43 

1075 

0.51 

3670 

PR  =  2.5 

0.69 

2500 

0.76 

4110 

PR  =  7.5 

0.44 

3300 

0.51 

8650 
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T1  T2  T3  T4  T5  T6  T7  T8  T9  T10  Til  T12  T13  T14  T15 
State  Points 

Fig.  8.  Temperature  distribution  of  oxygen  ion  and  hydrogen  proton-conducting  DA- 
SOFC/GT  system  (base  case). 


The  first  and  second  law  efficiencies  of  the  base  case  study  are 
listed  in  Table  5.  The  data  revealed  the  superiority  of  the  system 
with  SOFC-H  option  mainly  due  to  its  higher  exergy  efficiency 
when  compared  to  that  of  SOFC-O.  This  difference  is  expected  since 
SOFC-H  boasted  lower  ohmic  and  concentration  polarizations. 

5.3.1.  Effect  of  ammonia  utilization 

Fig.  10  reflects  the  effect  of  increasing  the  ammonia  utilization  in 
the  fuel  cell  module  on  the  total  specific  power  of  the  plant  and  the 
back  work  ratio  of  the  gas  turbine-compressor  assembly.  Although 
the  system  with  SOFC-H  benefited  from  higher  specific  power 
in  comparison  to  that  with  SOFC-O,  both  systems  experienced 
considerable  drop  in  specific  power  as  the  fuel  utilization  was 
increased.  Consequently,  the  decrease  in  the  molar  flow  rate  of  the 
gaseous  mixture  entering  the  gas  turbine  resulted  in  reducing 
the  production  of  mechanical  power.  However,  the  required 
compressor  power  was  also  reduced  at  a  different  rate  leading  to 
a  slight  increase  in  the  back  work  ratio. 

In  addition,  Fig.  11  emphasizes  the  decreasing  trend  of  energy 
and  exergy  efficiencies  of  the  integrated  system  in  response  to  the 
increase  of  ammonia  utilization.  It  is  therefore  very  important  to 
optimize  the  fuel  utilization  of  the  system  against  the  competing 
objective  of  maximum  specific  power. 

5.3.2.  Effect  of  fuel  cell  temperature 

Increasing  the  operating  temperature  of  the  fuel  cell  resulted  in 
appreciable  decline  of  the  total  specific  power  as  revealed  in  Fig.  12. 
At  elevated  temperatures,  the  system  achieved  lower  operating 
power  densities  to  allow  for  sufficient  heat  generation  to  satisfy  the 
heating  requirements  of  the  reactants  entering  the  fuel  cell. 
Conversely,  the  higher  temperature  of  the  gaseous  products 
increased  the  power  production  of  the  gas  turbine  hence  the  lower 
back  work  ratio. 


Table  4 

Specific  power  and  exergy  destruction  in  kj  mol1  (base  case). 


SOFC-O 

SOFC-H 

Overall  system 

Total  specific  power 

203.6 

219.5 

Total  exergy  destruction 

100.0 

83.5 

Component 

Specific  power  of  fuel  cell 

Specific  power  of  electric  generator 

126.0 

77.6 

141.9 

77.6 

65.56% 


^HX2 

^^HX3 

BMC 


17.66% 


7.69% 


6.89% 


0  41%  0.12%  168% 


E 


Fig.  9.  The  distribution  of  exergy  destruction  by  system  components  (base  case).  Part  A 
for  SOFC-O  option.  Part  B  for  SOFC-H  option. 


The  effect  of  the  fuel  cell  temperature  on  the  efficiency  of  the 
systems  is  shown  in  Fig.  13.  It  can  be  observed  that  the  system  with 
SOFC-H  option  reflected  superior  performance  and  efficiency  in 
comparison  to  that  coupled  with  SOFC-O.  As  the  temperature  of  the 
fuel  cell  increased,  the  electric  energy  and  exergy  efficiencies 
degraded  due  to  lower  specific  power.  Notwithstanding,  the  higher 
temperature  of  the  gaseous  products  boosted  the  amount  of  recov¬ 
erable  heat  and  the  associated  exergy  leading  to  an  overall  increase 
in  the  total  energy  and  exergy  efficiencies  of  the  respective  system. 


Table  5 

Energy  and  exergy  efficiencies  of  the  DA-SOFC/GT  system  (base  case). 


SOFC-O 

SOFC-H 

Overall  system 

V  t 

76.7 

81.1 

it 

69.9 

74.3 

Ve 

64.3 

69.3 

ie 

60.2 

65.0 

Component 

iFC 

66.7 

76.6 

icc 

93.6 

96.0 

ir 

95.1 

95.1 

ic 

90.2 

90.2 

iHx  1 

61.2 

61.2 

inxi 

96.4 

96.4 

iHx  3 

99.3 

99.3 

Efficiency  [%]  Efficiency  [%]  Specific  Power  [kj  mof’] 
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Ammonia  Utilization  [%] 


Temperarure  [K] 


l  Effect  of  ammonia  utilization  on  the  total  specific  power  and  back  work  ratio  of 
stem. 


Fig.  12.  The  effect  of  increasing  the  fuel  cell  temperature  on  the  total  specific  power 
and  back  work  ratio. 


Ammonia  Utilization  |%]  Temperature  (K) 


Fig.  11.  Effect  of  ammonia  utilization  on  the  energy  and  exergy  efficiency  of  the 
system.  Part  A  for  SOFC-O  option.  Part  B  for  SOFC-H  option. 


Fig.  13.  The  effect  of  increasing  the  fuel  cell  temperature  on  the  efficiency  of  the 
system.  Part  A  for  SOFC-O  option.  Part  B  for  SOFC-H  option. 
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Fig.  14.  The  effect  of  pressure  ratio  on  the  total  specific  power  and  back  work  ratio  of 
the  integrated  system. 


Fig.  15.  The  effect  of  the  pressure  ratio  on  the  efficiency  of  the  integrated  system.  Part 
A  for  SOFC-O  option.  Part  B  for  SOFC-H  option. 


Table  6 

The  effect  of  pressure  ratio  on  the  cooling  performance  of  ammonia  of  the  integrated 
system. 


Pressure  ratio 

2.5 

5.0 

7.5 

Tank  temperature  [K] 

259.5 

277.3 

289 

Cooling  load  [kj  mol-1] 

47.3 

46.7 

46.3 

Effectiveness  [%] 

14.9 

14.7 

14.6 

5.3.3.  Effect  of  pressure  ratio 

In  our  previous  work  [2],  we  demonstrated  that  increasing  the 
operating  pressure  enhanced  the  performance  of  solid  oxide  fuel 
cells.  However,  Fig.  14  reflects  the  adverse  impact  of  increasing  the 
operating  pressure  of  the  integrated  system  on  the  total  specific 
power  and  the  back  work  ratio.  This  indicates  that  a  significant 
portion  of  the  mechanical  power  produced  is  drawn  by  the  air 
compressor  thereby  reducing  the  net  electric  power  generated  per 
unit  mole  of  fuel.  Accordingly,  the  system  suffered  a  considerable 
decline  in  the  electric  energy  and  exergy  efficiencies  within  the 
studied  range  as  shown  in  Fig.  15.  Nonetheless,  the  ratio  of  recov¬ 
erable  heat  from  the  recuperator  to  that  added  to  the  ammonia  tank 
increased  proportionally  with  the  increase  of  the  pressure  ratio.  In 
this  manner,  the  reduction  of  the  overall  energy  and  exergy  effi¬ 
ciencies  was  less  pronounced  at  higher  pressure  ratios. 

It  is  also  worth  mentioning  that  the  cooling  performance  of  the 
system  is  strongly  affected  by  the  pressure  ratio  as  demonstrated  in 
Table  6.  Although  the  cooling  load  capacity  and  effectiveness 
remained  virtually  unchanged,  the  refrigeration  or  cooling  poten¬ 
tial  may  be  diminished  due  to  negative  impact  of  increasing  the 
operating  pressure  of  the  system  on  the  saturation  temperature  of 
ammonia. 

5.3.4.  Effect  of  reference  state 

Exergy  is  a  property  of  the  system  and  its  environment.  Fig.  16 
depicts  the  effect  of  changing  the  reference  temperature  on  the 
exergy  efficiencies  of  the  integrated  system.  While  the  effect  is 
negligible  on  the  exergy  efficiency  based  on  electric  power 
production,  the  total  exergy  efficiency  suffers  a  minor  drop  over 
the  studied  range.  The  reduced  efficiency  is  proportional  to  the 
decrease  in  exergy  associated  with  the  heat  recovered  from  the 
recuperator.  This  exergy  is  governed  by  a  factor  corresponds  to 


Fig.  16.  Effect  of  reference  temperature  on  the  exergy  efficiency  of  the  system  (base 
case). 


84 


F.  Ishak  et  al.  /  Journal  of  Power  Sources  212  (2012)  73-85 


efficiency  of  a  Carnot  engine  operating  between  the  temperature  of 
the  gaseous  stream  and  the  reference  temperature.  Hence, 
increasing  the  reference  temperature  reduces  the  amount  of  useful 
work  potential  of  the  heat  recovery  process  resulting  in  lower  total 
exergy  efficiency. 

6.  Conclusions 

This  study  has  presented  an  innovative  hybrid  system, 
comprising  of  ammonia-fed  solid  oxide  fuel  cell  integrated  with 
a  gas  turbine  power  cycle.  The  simple  yet  effective  configuration  of 
the  system  allows  makes  it  advantageous  in  combined  heating, 
cooling  and  power  applications. 

The  thermodynamic  and  the  electrochemical  processes  of  the 
system  are  modelled  and  validated,  with  excellent  agreement, 
against  experimental  data  obtained  from  the  literature.  In  addition, 
a  parametric  investigation  is  performed  to  examine  the  effects  of 
design  and  operating  variables  on  the  performance  and  thermo¬ 
dynamic  efficiency  of  the  system.  It  is  evident  from  the  results  of 
the  study  that  the  overall  system  performance  is  affected  by  the 
configuration  and  arrangement  of  system  components  among 
other  factors. 

The  highest  exergy  destruction  found  in  the  fuel  cell  can  be 
attributed  to  the  thermodynamic  and  electrochemical  irrevers¬ 
ibilities.  The  second  highest  exergy  destruction  occurs  in  heat 
exchanger  (HX1),  mainly  due  temperature  mismatch  between  the 
incoming  fluid  streams. 

In  summary  the  hydrogen  proton-conducting  DA-SOFC 
demonstrates  a  better  performance  in  comparison  to  the  ion¬ 
conducting  counterpart  due  to  higher  partial  pressure  of 
hydrogen  at  the  anode  under  all  operating  conditions.  Nonetheless, 
it  is  imperative  to  optimize  all  design  and  operating  parameters  to 
achieve  the  best  possible  performance  and  efficiency  of  the  system. 

Nomenclature 

A  Area  [m2] 

Bo  Permeability  [m2] 

cp  Specific  molar  heat  capacity  [J  mol-1  I<-1] 

dx  Depth  of  electrode  layer  [m] 

Deff  Effective  binary  diffusion  [m2  s-1] 

D£ff  Effective  Knudsen  diffusion  [m2  s-1] 
ex  Molar  specific  exergy  [J  mol-1] 

E  Reversible  cell  potential  or  voltage  [V] 

E°  Standard  cell  potential  or  voltage  [V] 

Ex  Exergy  rate  [W] 

Exd  Exergy  destruction  rate  [W] 

f  Fugacity  of  chemical  species  in  the  system 

f°  Fugacity  of  chemical  species  at  STP 

F  Faraday’s  constant  (96,485  C  mol-1) 

g  Specific  molar  Gibbs  energy  [J  mol-1] 

G  Gibbs  energy  [J] 

g°  Specific  molar  Gibbs  energy  at  STP  [J  mol-1] 

h  Specific  molar  enthalpy  [J  mol-1] 

H  Enthalpy  [J] 

H°  Standard  Enthalpy  at  STP  [J] 

Hc  Heat  capacity  rate  [W  I<-1] 

I  Electric  current  [A] 

J  Current  density  [A  m-2] 

Jo  Exchange  current  density  [A  m-2] 

k  Ratio  of  specific  molar  heat 

M  Molar  mass  [g  mol-1] 

n  Number  of  moles  [mol] 

h  Molar  flow  rate  [mol  s-1] 

N  Molar  flux  [mol  m-2  s-1] 

p  Partial  pressure  [kPa] 


p 

Total  pressure  [kPa] 

PR 

Pressure  ratio 

Q 

Heat  [j] 

Q 

Rate  of  heat  transfer  [W] 

rP 

Mean  pore  radius  of  electrode  [m] 

R 

Gas  constant  [J  mol-1  K-1] 

Rq 

Electrolyte  resistivity  [Q  m] 

s 

Molar  specific  entropy  [J  mol-1  K-1] 

S 

Entropy  [J  K-1] 

T 

Temperature  [K] 

U 

Utilization  [%] 

V 

Specific  molar  volume  [m3  mol-1] 

V 

Working  potential  or  voltage  [V] 

W 

Work  [J] 

W 

Work  rate,  power  [W] 

X 

Molar  ratio  of  produced  hydrogen  to  consumed  fuel 

y 

Molar  fraction  of  chemical  species 

z 

Number  of  electrons  transferred  per  mole  of  fuel 

Greek  Letters 

a 

Charge  transfer  coefficient 

r 

Dimensionless  temperature 

5 

Electrolyte  thickness  [m] 

A 

Net  change  of  quantity 

e 

Porosity  of  electrode 

e\x 

Lennard-Jones  temperature  parameter  [I<] 

V 

Energy  or  thermal  efficiency 

P 

Viscosity  [N  s  m-2] 

? 

Heat  exchanger  effectiveness 

a 

Molecular  collision  diameter  [A] 

T 

Tortuosity  of  electrode 

(Pact,  an 

Anode  activation  over-potential  [V] 

(Pact,  ca 

Cathode  activation  over-potential  [V] 

(Pconc,  an 

Anode  concentration  over-potential  [V] 

(Pconc,  ca 

Cathode  concentration  over-potential  [V] 

Ohmic  over-potential  of  electrolyte  [V] 

* 

Exergy  efficiency 

Qd 

Collision  integral  of  diffusion 

Subscripts 

AC 

Alternating  current 

amb 

Ambient 

c 

Cold 

ch 

Chemical 

C 

Compressor 

CC 

Combustion  chamber 

CL 

Cooling  load 

cmp 

Component 

DC 

Direct  current 

e 

Electric 

eff 

Effective 

/ 

Fuel 

FC 

Fuel  cell 

G 

Generator 

h 

Hot 

HX 

Heat  exchanger 

i 

Chemical  species 

in 

Inlet 

INV 

Inverter 

j 

Chemical  species,  system  component 

MC 

Mixing  chamber 

min 

Minimum 

mix 

Mixture 

k 

Knudsen  diffusion 

0 

Oxidant 
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out  Outlet 

T  Turbine 

ph  physical 

r  Reactants 

REC  Recovery,  Recuperator 

s  Isentropic 

T  Turbine 

Acronyms 

BWR  Back  Work  Ratio 

CHP  Combined  Heat  and  Power 

CHCP  Combined  Heating,  Cooling  and  Power 

DA-SOFC  Direct  Ammonia  Solid  Oxide  Fuel  Cell 
DA-SOFC/GT  Direct  Ammonia  Solid  Oxide  Fuel  Cell  Integrated  with 
Gas  Turbine 
DGM  Dusty  Gas  Model 

GEM  Gibbs  Energy  Minimization 

LHV  Lower  Heating  Value 

LSM-YSZ  Lanthanum  Strontium  Manganite/Yttria-stabilized 
Zirconia 

Ni-YSZ  Nickel/Yttria-stabilized  Zirconia 

SOFC/GT  Solid  Oxide  Fuel  Cell  Integrated  with  Gas  Turbine 
SOFC-H  Hydrogen  proton  (H+)  conducting  Solid  Oxide  Fuel  Cell 
SOFC-O  Oxygen  ion  (02~)  conducting  Solid  Oxide  Fuel  Cell 
STP  Standard  Temperature  and  Pressure,  298  I<  and 

101.325  kPa  respectively 
TPB  Triple-Phase  Boundary 

YSZ  Yttria-stabilized  Zirconia 
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